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Overlap and domain wall fermions: what is the price of chirality? 
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In this contribution the costs of simulations employing domain wall and overlap fermions are estimated. In the 
discussion we will stay within the quenched approximation. 



Chirally invariant formulations of lattice QCD 
are still relatively new. However, besides their 
conceptual advantages, a number of physics appli- 
cations have already been performed jlj] . In these 
works it became clear that simulations with such 
formulations are expensive, even when they are 
restricted to the quenched approximation. Con- 
sequently, at the time of the conference no simu- 
lations have been performed with such dynamical 
quarks so far. 

Our conclusions concerning the costs of sim- 
ulations with dynamical overlap or domain wall 
fermions will therefore be rather indirect: we will 
estimate the overhead of using chiral invariant 
formulations to standard Wilson fermions in the 
quenched case and assume this overhead to be the 
same for an unquenched simulation. The discus- 
sion below will be organized in the form of four 
statements. 

Statement 1 

The 5- dimensional domain wall construction 
in the limit N s — > oo is completely equivalent 
to a ^-dimensional lattice formulation of overlap 
fermions. 

This is a mathematical statement that can be 
proven rigorously Denoting by s the extra 

5th dimension, the 5-dimensional operator reads 



D 5 = 



d s ) - a s d* 8 d s ] + M 



with 

M = D,„ 



to - 



(1) 



(2) 



Here D w is the standard Wilson-Dirac operator, 
a s the lattice spacing in the extra dimension and 
Too a mass parameter. 

Keeping a s finite, in the limit N s — > oo we ob- 
tain a 4-dimensional operator 



aDi = 1 - A(A^A)- 1 / 2 
with 

A = -a s M(2 + a s My 1 



(3) 



(4) 



that satisfies the Ginsparg- Wilson relation. If 
also a s — > 0, Neuberger's overlap operator j| is 
recovered. 

Statement 2 

Domain wall fermions do not perform better 
than overlap fermions. 

The mathematical equivalence of domain wall 
and overlap fermions gives rise to the suspicion 
that also in practical applications no clear prefer- 
ence for either formulation as far as the expense 
of the simulation is concerned- can be given. We 
give one example below for the cost of comput- 
ing the pion propagator rV to a certain relative 
precision. In fig. 1 we plot the ratio R, 



R 



fc (T/3) -T a PP rox (T/3) 



r° xact (T/3) 



(5) 
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Aiming at a, say, per mil precision we find for 
this particular example that domain wall fermions 
are about a factor of two more expensive than 
overlap fermions. It is clear that for both for- 
mulations additional improvements might be im- 
plemented. However, it seems very unlikely that 
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one particular formulation will give an order of 
magnitude better performance. 



10 
1 

0.1 



0.01 r 

0.001 r 

0.0001 r 
lO- 5 

io- 6 b 
io- 7 



■ GW 

• DW 



100 200 300 400 
time 



Figure 1 . The ration R of eq. (|^) as a function of 
time needed (on a particular machine) to obtain a 
relativ precision for the pion propagator. Results 
are averaged over ten configurations on a 8 3 x 24 
lattice at /? = 5.85. 

Statement 3 

Keeping N s finite: the residual mass is not all. 

The locality of chirality breaking effects is only 
guaranteed at distances \x — y\ 3> N s a. This 
is discussed in First investigations of this 

question have been performed in || . I find it very 
important to study the effects of this observation 
further on a quantitative level. 

Statement 4 

Whatever time estimate is found for Wilson 
fermions: multiply the effort by a factor of 
O(100) for chiral symmetric actions. 

The factor referred to in this statement is really 
only an order of magnitude estimate. The reasons 
for this large factor are different for overlap and 
domain wall fermions. 

In the case of overlap fermions it turns out 
that the polynomial required to approximate the 
square root has typically a degree of O(100). 



Since this polynomial has to be evaluated in every 
step of a linear solver, the cost of simulations with 
overlap fermions increase correspondingly. It is 
also noteworthy that so far no way of precondi- 
tioning the overlap operator has been found. 

In the case of domain wall fermions there is 
some experience that the number of conjugate 
gradient iterations is larger for the 5-dimensional 
problem. This in addition to the additional num- 
ber of slices in the extra dimension gives again 
a large value of the factor compared to standard 
Wilson fermions. 

Conclusion 

If the estimate of about lOOTflops to solve most 
of the problems in lattice QCD as discussed in 
the panel contributions is indeed correct, then 
this would mean a demand of 10 Petaflops for 
simulations with chirally invariant formulations 
of lattice-QCD. We then would be close to the 
number |hJ anticipated by K. Wilson in his 1989 
Capri contribution. 

When, as it is done at this conference, we em- 
phasize that thinking about better algorithms is 
one of the most important things the lattice com- 
munity should address, then this is even more 
true for chirally invariant formulations of lattice- 
QCD. And, finally, the hope is that putting effort 
into the development of better algorithms, it can 
be shown that statement number 4 is not correct. 
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